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Hydrogen adatoms are one of the most the promising proposals for the functionalization of 
graphene. Hydrogen induces narrow resonances near the Dirac energy, which lead to the forma¬ 
tion of magnetic moments. Furthermore, they also create local lattice distortions which enhance the 
spin-orbit coupling. The combination of magnetism and spin-orbit coupling allows for a rich variety 
of phases, some of which have non trivial topological features. We analyze the interplay between 
magnetism and spin-orbit coupling in ordered arrays of hydrogen on graphene monolayers, and clas¬ 
sify the different phases that may arise. We extend our model to consider arrays of adsorbates in 
graphene-like crystals with stronger intrinsic spin-orbit couplings. 

PACS numbers: 


I. INTRODUCTION 

Adsorbates and vacancy effects in graphene have been 
a major field of research in recent years. An isolated va¬ 
cancy gives rise to localized resonance states near the 
Fermi level [U [2]. Hydrogen adatoms are expected to 
form a strong covalent bond with a carbon atom in the 
graphene lattice. The bond effectively removes one tt or¬ 
bital in the graphene band, leading to a sharp resonance 
near the Fermi energy, in a similar way to the case of 
the vacancy [3]. The electron-electron repulsion prevents 
doubly occupied states of this resonance, and lead to the 
formation of a magnetic moment. In addition, the carbon 
atom coupled to the hydrogen atom is displaced from the 
graphene plane, inducing a local sp^ hybridization, which 
increases the spin-orbit coupling [H |5]. 

The intrinsic spin-orbit coupling in a perfect graphene 
layer creates a gap, and turns graphene into a topological 
insulator [6]. The presence of magnetic moments induces 
a exchange coupling with the spins of itinerant electrons, 
and breaks time reversal symmetry. The combination of 
a uniform exchange coupling and the extrinsic Rashba 
coupling in graphene leads to a quantum anomalous Hall 
phase [7]. This phase is an example of systems which do 
not show time reversal symmetry and have topologically 
protected edge states without Landau levels [8]. 

The effect of a uniform magnetic field and the spin- 
orbit coupling has been extensively studied in silicene[9]- 
m- The possibility of inducing non trivial topological 
features in the electronic structure of graphene by the 
addition of heavy atoms has also been considered [12]- 
[14], and non trivial features associated to the spin- 
orbit coupling have been found in lead intercalated GVD 
graphene [21]. 

On one hand, the formation of magnetic moments near 
vacancies and hydrogen atoms in graphene has been ex¬ 
tensively investigated [15], both theoretically (see, for in¬ 
stance [T6ll2Q] ). and experimentally [221126] . On the other 
hand, it is known that adsorbates on graphene form or¬ 


dered arrays in a variety of situations [27l 1^ . Adatoms 
on graphene interact among themselves, and can form a 
variety of ordered patterns [291131] . More importantly, as 
shown recently, hydrogen adsorbates can be manipulated 
with an STM tip [32]. 

We analyze here the combined effects of the exchange 
coupling and modified spin-orbit coupling due to hydro¬ 
gen adsorbates which form a regular array. As men¬ 
tioned before, the formation of ordered arrays can be 
induced by the interactions between adatoms, or they 
can be fabricated on purpose. The results are expected 
to give hints on the local changes in the electronic struc¬ 
ture near an isolated impurity. In section II we present 
the model, and discuss the main features of arrays of 
hydrogen adatoms. In section III we study the topolog¬ 
ical properties of the hydrogenated graphene superlat¬ 
tice considering only the local adatom-induced SO cou¬ 
plings. In section IV we proceed to study related situa¬ 
tions, where the spin-orbit coupling throughout the en¬ 
tire lattice cannot be neglected. Finally, we discuss the 
most relevant results, and the open questions raised by 
our work. 


II. THE MODEL. 

We will study the topological properties of a periodic 
array composed of a hydrogen adatom on top of a carbon 
atom in the 5x5 supercell of graphene as shown in the 
first two panels of|n]. We describe the electronic struc¬ 
ture with a tight binding model with one Pz orbital per 
carbon atom and a nearest neighbor hopping parameter, 
t. The effect of the covalent bonding between hydrogen 
and carbon is approximated by a large shift of the energy 
of the Pz orbital nearest to the hydrogen atom, eo- For 
eo ^ \t\ a sharp resonance appears near the vacancy [T], 
which decays slowly as function of the distance. The res¬ 
onance is mostly localized in the sublattice which does 
not include the perturbed Pz orbital. 
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We assume that the electron-electron interaction in¬ 
duces a exchange splitting, between the spin up and 
spin down levels at the sites where the resonance is lo¬ 
cated. This is the effect which arises if the interaction 
is local and it is treated within a mean field approxima¬ 
tion. The orientation of the exchange field with respect 
to the plane of the graphene layer is determined by the 
spin-orbit coupling. To a first approximation, a Rashba- 
like coupling does not modify the spin, as a Rashba cou¬ 
pling implies hopping between different sublattices and 
the spin is localized in one sublattice (note that the itin¬ 
erant electrons are affected). The full Hamiltonian for 
the adatom in the super cell is: 


H = Hkin + Hz + Hso ( 1 ) 



where Hkin represents the kinetic terms of the Hamilto¬ 
nian, Hso represents the spin-orbit coupling terms, and 
Hz is the exchange splitting due to the magnetism in¬ 
duced by the onsite potential of the adatom, which takes 
the form of a usual site dependent Zeeman term. In this 
section we will explain the kinetic and exchange contri¬ 
butions to the Hamiltonian. In the next two sections we 
will discuss two cases separately: local SO coupling in¬ 
duced by the adatom, and uniform SO coupling induced 
by a substrate. 

The kinetic term can be written: 


Hkin = X] + ^o^cl;„^aCc„,s (2) 

{hj),s s 


where t is the hopping between nearest neighbors i and 
j in the honeycomb lattice, eo is a large onsite potential 
in the hydrogenated carbon site (in our calculation we 
consider it to be 100 times the hopping integral), and 
^Ch s {^Ch,s) creates (annihilates) an electron in the hy¬ 
drogenated carbon site {Ch)- The onsite potential due 
to the adatom gives rise to a gap between the conduction 
and valence bands, as well as a vacancy state lying at 
the Fermi energy. The size of this gap is proportional to 
the impurity concentration, so that in the dilute limit of 
an isolated vacancy in graphene, there is a zero energy 
vacancy state, but no gap. 

The vacancy state also produces a non-uniform in¬ 
ternal magnetization in the sub lattice opposite to the 
adatom, which decreases with distance from the impurity 
site. The local magnetization splits the impurity bands, 
creating a gap at the Fermi level, which can be seen in 
panel (d) of[H| We model this magnetization effect by 
the exchange splitting term in the Hamiltonian, which 
has the form: 


FIG. 1: Single hydrogen adatom in graphene sheet (a), and 
unit cell of the simulation (b). As shown in (b), Rashba and 
exchange fields are located only around the adatom. Band 
structure of a pristine 5x5 graphene unit cell (c), and with a 
vacancy and local exchange field (d). 





FIG. 2: Band structures (a,b) and Berry curvatures (c,d) 
of a graphene supercell with one hydrogen atom, with lo¬ 
cal Rashba field and local off-plane exchange. Panels (a,c) 
corresponds to realistic values of Xr, whereas panels (b,d) 
corresponds to a large Xr. 


potential in the B sites. The splitting magnitude is 
proportional to the magnitude of the vacancy state at 
site is, that is: 


Hz = ( 3 ) 

is 


AiB = (4) 


where the is label means that the summation is carried 
out over all sites of the B sub lattice opposite to Ch- 
Effectively, this term contributes a spin-dependent onsite 


After solving for the band structure of the spin- 
unpolarized kinetic part of the Hamiltonian, we introduce 
the sub lattice magnetization via eq. (©■ 
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III. HYDROGENATED GRAPHENE: LOCAL 
SPIN-ORBIT COUPLINGS 


In this section we discuss the model with local, 
adatom-induced SO couplings. The Hamiltonian for the 
spin orbit coupling induced by the adatom is [5]: 


Hso = 


+ 

+ 


( 5 ) 


3 E 

{{Ch,})) 

2i 
~3 


LV3 




+ H.C. 


- W c ■ . 


{{ChJ}} 


2i 


^ A 




2 / V 

(dd)) 


Xr{s X dcH,j)z 
X^PIa{s X Dij)z 


^H.c. 


Where dcnj is the unit vector connecting Ch to its near¬ 
est neighbors, and Dij is the unit vector connecting the 
second neighbors of Ch- The operator <2^^ ^ ^ ( Cnn,j,s) 
creates (annihilates) electrons of spin s in the jth near¬ 
est neighbor of Ch^ whereas c]j^ ^ creates an electron at 
(7i:f.The first term describes the adatom modified intrin¬ 
sic SO coupling A/, which induces spin-preserving hop¬ 
pings between Ch and its second neighbors. The sec¬ 
ond term {Xr) describes the Bychkov-Rashba hoppings 
caused by the local breaking of space inversion symmetry 
around the adatom site. It induces spin-flipping hoppings 
between Ch and its nearest neighbors. The third term 
{X^pia) induces spin-flipping hoppings between the near¬ 
est neighbors of Ch- We have neglected in this section 
the intrinsic SO coupling of pristine graphene, which in¬ 
duces hoppings between second neighbors not containing 
Ch- 

The presence of Rashba interaction and an exchange 
field in the sites where the impurity bands are located, 
strongly suggests the possibility of reaching a Quantum 
Anomalous Hall (QAH) state, via a topological phase 
transition occurring in the corresponding gap. We in¬ 
spect the topological properties of the valence bands of 
the bulk system, by calculating the Chern number 




^n{}^x-)^y)d k 


(6) 


where the integral is carried out over the first Brillouin 
zone, and the summation is carried out over all valence 
bands. The Chern number gives the number of chi¬ 
ral edge states at each edge of a nano ribbon, and is 
related to the quantized charge Hall conductance via 
= Ce^Ih . Using the ab-initio values reported by 
[5] for the spin-orbit coupling strengths and onsite po¬ 
tential, we calculate the Berry curvature of the graphene 
supercell with an adatom. The SO coupling strengths 
are Xr = 1.14 x 10“^ t, X^ria = —2.66 x 10“^ t and 
A/ = —7.26 X 10“^ t. The magnitude of the exchange 
splitting at the nearest neighbors of the adatom site is 
set at 0.1 t throughout this calculation. A close-up of 


the band structure around the gap, and the Berry curva¬ 
ture profiles for these values are shown in panels (a) and 

(c) of Fig. (2), integration over the Brillouin zone yields 
C = 0. The Berry curvature has different signs in op¬ 
posite valleys and the Chern number around each valley 
is not an integer. This situation is similar to the Valley 
Hall effect in graphene in the presence of a sub lattice 
staggered potential, where valley currents are formed in 
the direction transverse to an applied in-plane electric 
field. 

Moving away from graphene, we will explore this model 
in the regime of stronger SOC, which will be repre¬ 
sentative of heavier graphene like honeycomb lattices. 
We will investigate the conditions for achieving a QAH 
state in this setting, we do this by checking the low- 
energy band structure and Berry curvature while increas¬ 
ing the Rashba parameter. A closing and reopening of 
the gap occurs at the topological phase transition from 
the Valley-Hall regime to the QAH state. We found this 
to occur for values of the Rashba parameter that are 10^ 
times larger than the ab-initio values of [5] {Xr = 3.25 t). 
Panel (b) in Fig. (2) shows the band structure for this 
situation of very large Rashba coupling, whereas panel 

(d) shows the profile of the Berry curvature in the same 
situation. The magnetization breaks both time reversal 
symmetry and sub lattice symmetry, since the adatoms 
are located only in one sublattice. We find the Berry 
curvature to be almost entirely concentrated around the 
K valley. We find C = 1. We have investigated the band 
gaps at half-filling for lower adatom concentrations, and 
we have found the QAH phase to exist within a region 
of parameter space beyond physically realistic values [33]. 
Therefore, although theoretically the simultaneous ex¬ 
istence of exchange and Rashba couplings in graphene 
might be able to create a quantum anomalous Hall state, 
within the realistic values such phase is not expected 
to be observed. Nevertheless, it is worth to note that 
the present model can be also applied to model further 
2D honeycomb systems with stronger SOC, in which the 
topological state could be realized. In the following sec¬ 
tion, we will show how in these systems, where apart form 
the Rashba coupling there is a strong intrinsic SOC, the 
topological state can be easily realized. 


IV. HONEYCOMB LATTICE WITH 
ADSORBATES AND UNIFORM SOC 


We have seen that the hydrogen-induced local Rashba 
is not strong enough to induce a topological phase tran¬ 
sition to a QAH insulator, now we investigate the possi¬ 
bility of achieving such a state in a different setting: a 
supercell of a graphene like honeycomb crystal in which 
the intrinsic SO coupling cannot be neglected, in the 
presence of a substrate that induces a uniform Rashba 
effect, see panel (a) in Fig. (3). This model can be 
applied to a large family of systems such as silicene, ger- 
manene, stanene, hydrogenated bismuth, metal organic 
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FIG. 3: Scheme (a) of unit cell with vacancy induced local 
magnetism, and uniform Rashba induced by off-plane electric 
field. Band structure with vacancy, local magnetism, uniform 
intrinsic SOC with (b) and without Rashba (c). Berry curva¬ 
ture for (b), showing a Chern number C — -\-l. In comparison 
with Fig.2, the anomalous Hall state (b) appears at arbitrary 
small Rashba field, provided intrinsic SOC closes the band 
gap (c). 


frameworks [35H4Q] .In this model, the role of the adsor¬ 
bate impurity is to generate the fiat impurity band and 
to induce a non-uniform exchange field. The spin-orbit 
coupling Hamiltonian in this model has the form: 

Hso = (7) 

2i + 

+ 

(hi) 

+ I 

where the SO hoppings are now uniformly distributed 
around the supercell. Apart from this, the supercell 
Hamiltonian remains unchanged. Since the magnitude 
of the gap between the impurity bands decreases with 
adatom concentration, we can consider a regime in which 
the Kane-Mele intrinsic SO coupling is strong enough to 
close the exchange-splitted gap, as shown in panel (c) of 
Fig. (3). This situation can be observed with A/ = 0.02 
t, Xr = 0.0 t, and a exchange splitting of 0.25 t in the 
nearest-neighbors of the adatom. There is no gap around 
the K valley in this regime. Once the Rashba SO coupling 
is induced, a gap with topologically non-trivial proper¬ 
ties appears as shown in panels (b) and (d) of Fig. (3), 
for the band structure and Berry curvature respectively, 
for a Rashba of Xr = 0.02 t. Integration yields a Chern 
number C = 1, signaling a QAH state with 1 edge state. 

A qualitative description of the properties of the sys¬ 
tem can be seen from the phase diagram in panel (a) of 
Fig. (4). For small values of the intrinsic and Rashba SO 


couplings the system is a trivial insulator. By increasing 
the values of either the Kane-Mele or Rashba the system 
turns into a metal, whereas the combination of higher 
intrinsic and Rashba couplings changes the system into 
a QAH insulator with one conducting edge state. The 
regime where the system is a trivial insulator, is charac¬ 
terized by Berry curvatures of opposite signs in opposite 
valleys. 

So far, we have considered the topological properties of 
the superlattice of adsorbs at half-filling. The adsorbates 
break sub lattice inversion symmetry, and in the trivial 
insulating case, they generate a Berry curvature distri¬ 
bution similar to that of graphene in the presence of a 
sublattice staggered potential, characterized by peaks of 
opposite signs at the two valleys. When the Fermi level is 
tuned to be at the bottom of the conduction band, which 
can be done by applying a gate voltage, the Berry curva¬ 
ture distribution continues to be asymmetric in opposite 
valleys, as shown in panels (c) and (d) of Fig. (4)., giv¬ 
ing rise to a non-quantized Valley-Hall effect, i.e. skew¬ 
scattering for electrons in opposite valleys (see panel b), 
and thus a net valley current that travels in the bulk. The 
sign of the anomalous velocity can be switched by chang¬ 
ing the sub lattice in which the adsorbates are placed. 

This later effect, independent of the Rashba coupling, 
is the one to be observed in graphene, provided all the 
hydrogen impurities are located in the same sublattice. 
Nevertheless, in a realistic random hydrogenation, pro¬ 
vided clustering can be avoided im, in a simplistic pic¬ 
ture it is expected that the hydrogen atom will lie ran¬ 
domly in both sublattices, therefore giving rise to vanish¬ 
ing valley Hall effect. However, if all the hydrogen atoms 
are located in the same sublattice, the net valley Hall 
current will be different from zero and strong non-local 
signal might be observed, very much like in gated bilayer 
graphene [42] , 


V. SUMMARY AND DISCUSSION 

We have aimed at an understanding of the topological 
properties of hydrogenated graphene in the dilute limit 
by a Berry curvature analysis of the bulk gap at half 
filling. This gap is opened by the splitting of the im¬ 
purity bands due to the exchange field. It is influenced 
by adatom concentration inasmuch as the vacancy state 
and the induced magnetic moment depend on the dis¬ 
tance between adatoms [43]. We have found that within 
the range of realistic values for the SO couplings and on¬ 
site potential reported in [5] the sub lattice asymmetry 
induced by the adatom predominates over the exchange 
field and local SO couplings. The predominant phase 
is similar to the Valley Hall effect in gapped graphene, 
where electrons in opposite valleys acquire anomalous ve¬ 
locities transverse to an in-plane electric field in opposite 
directions, generating thus a valley current that travels 
in the bulk. We have also found that this system can 
theoretically undergo a topological phase transition to 




A/ 

V3 


H.c. 
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FIG. 4: (a) Schematic phase diagram of the adatom in the 
supercell with intrinsic and Rashba SO couplings as parame¬ 
ters, the exact values of the SO couplings vary as a function 
of impurity concentration, (b) The Berry curvature is con¬ 
centrated in one valley or the other depending on which sub 
lattice the impurities are located. Valley currents of opposite 
signs may be thus be engineered at will, (c) and (d) show re¬ 
spectively, the band structure and Berry curvature for slightly 
doped graphene. If the Fermi level now lies at the bottom of 
the conduction band the Berry curvature still has opposite 
signs at opposite valleys. 


a valley-polarized Quantum Anomalous Hall state, but 
this QAH state lives in a region of parameter space with 


Rashba coupling values that are 10^ times larger than 
the ab-initio values reported by [5]. 

We also extended our model in order to study the topo¬ 
logical properties of arrays of adsorbates in honeycomb 
crystals in which the uniform intrinsic and Rashba SO 
couplings cannot be neglected, such as stanene, silecene 
or germanene. Starting from a situation in which the 
intrinsic SOC is strong enough to mix the exchange- 
splitted impurity bands, a situation that can be realized 
for low concentrations of impurities, we have found that 
the Rashba parameter is capable of opening a gap with 
non-trivial topological properties, effectively turning the 
system into a QAH insulator. Finally we also consid¬ 
ered the situation in which the Fermi level lies inside the 
bottom of the conduction band, we have found a Berry 
curvature distribution of peaks of opposite signs in the 
two valleys, marking the the presence of bulk valley cur¬ 
rents in a direction transverse to an applied electric field. 
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